Diffusion and Current of Brownian Particles in Tilted Piecewise Linear Potentials: 

Amplification and Coherence 
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Overdamped motion of Brownian particles in tilted piecewise linear periodic potentials is con- 
sidered. Explicit algebraic expressions for the diffusion coefficient, current, and coherence level of 
Brownian transport are derived. Their dependencies on temperature, tilting force, and the shape of 
the potential are analyzed. The necessary and sufficient conditions for the non-monotonic behavior 
of the diffusion coefficient as a function of temperature are determined. The diffusion coefficient 
and coherence level are found to be extremely sensitive to the asymmetry of the potential. It is 
established that at the values of the external force, for which the enhancement of diffusion is most 
rapid, the level of coherence has a wide plateau at low temperatures with the value of the Peclet 
factor 2. An interpretation of the amplification of diffusion in comparison with free thermal diffusion 
in terms of probability distribution is proposed. 

PACS numbers: 05.40.-a, 05.60.-k, 02.50.Ey 



I. INTRODUCTION 



It is difficult to overestimate the importance of the 
effects caused by Brownian motion for soft condensed 
matter physics. An object of special attention has been 
Brownian motion in periodic structures which has var- 
ious applications in condensed matter physics, chem- 
ical physics, nanotechnology, and molecular biology 

Furthermore, the idea that adding noise to 
deterministic motion can give nontrivial results has led to 
many important discoveries, such as stochastic resonance 
[3 , resonant activation , noise-induced spatial patterns 
[31, noise- induced multistability as well as noise-induced 
phase transitions jS^ A A .10, .11.] , ratchets 2] , and hy- 
persensitive transport |l2| , to name but a few of the new 
phenomena in this field. 



The acceleration of thermal diffusion of Brownian par- 
ticles in tilted periodic potentials has been studied re- 
cently in a number of papers IS El El El. An smaU 
enhancement of diffusion induced by bias with respect 
to free diffusion was established for the asymptotic lim- 
its Is] . The papers Ti', T?, provided exact analytic 
treatments of the effect for arbitrary temperature, tilting 
force, and periodic potential. A giant amplification of dif- 
fusion up to fourteen orders of magnitude was predicted 
in Refs. 14. .ISj. Another new effect, a non-monotonic 
behavior of the diffusion coefficient and coherence level 
of the transport of Brownian particles as a function of 
temperature was found in Ref. Similar anomalies 

were observed in systems with spatially periodic temper- 
ature . Recently it was shown that non- homogeneous 
dissipation can induce enhancement and suppression of 
the diffusion as a function of temperature, as well as an 
increase of the coherence level of Brownian motion in 
tilted symmetric periodic potential Ell ■ A distinctive be- 
havior of diffusion also occurs if one applies a sinusoidal 
time-periodic force to the system, in addition to a tilted 
sinusoidal space-periodic force. It was demonstrated nu- 
merically that the interplay between frequency-locking 



and noise gives rise to a multi-enhancement of the effec- 
tive diffusion, and a rich behavior of the current, includ- 
ing partial suppression and characteristic resonances E3 
(see also 



Usually, when addressing diffusion enhancement, tilted 
harmonic spatial periodic potential is used, confer how- 
ever Ref. Eg- Furthermore, it is known that some fea- 
tures of the ratchet transport mechanism (e.g., the cur- 
rent reversals) are extremely sensitive to the shape of 
the potential j2l|. In the present paper, we carry out a 
comprehensive study of the dependence of diffusive and 
coherent motion of overdamped Brownian particles on 
temperature and tilting force for various shapes of tilted 
piecewise linear periodic potentials. The latter are suffi- 
ciently simple to allow an algebraic treatment of the rel- 
evant quantities, being at the same time physically rich 
enough to provide most of the effects characteristic of 
tilted periodic potentials with one minimum per period. 



The paper is organized as follows. Proceeding from 
the general scheme developed in Refs. 0, ^| (outlined 
in Sec. II), we derive in Sec. Ill the algebraic expressions 
for the diffusion coefficient, current, and Peclet factor 
for the case of a general tilted piecewise linear periodic 
potential. The asymptotic limits and particular cases are 
considered analytically in Sec. IV. In Sec. V we analyze 
the behavior of the diffusion coefficient and Peclet factor 
in the space of system parameters, i. e., the tilting force, 
temperature, and asymmetry parameter of the potential. 
We also discuss the effect of the amplification of diffusion 
in terms of probability distribution. The existence of 
correlation between the enhancement of diffusion and the 
stabilization of the coherence level of Brownian transport 
is demonstrated in Sec. VI. Our results are summarized 
in Sec. VII. 
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II. GENERAL SCHEME 



III. ANALYTIC COMPUTATIONS 



We consider overdamped motion of Brownian particles 
described by one-dimensional Langevin equation 



V- 



dx{t) _ dV{x) 



dt 



dx 



(1) 



where V{x) = Vq{x) - Fx and Vo{x) = Vo{x + L) is the 
periodic potential, F is the static tilting force, 77 is the 
coefficient of viscous friction, £^{t) is the Gaussian white 
noise with the mean value (^(i))=0 and the correlation 
function {^{t)^{t'))=2T]kBTS{t^t'). 

As usual, the effective diffusion coefficient is defined by 



D = lim 



2t 



In the same time-scale the current is determined as 
(i). lim Ml. 

t — >oo t 



(2) 



(3) 



According to Refs. the diffusion coefficient for the 

model O is equal to {F > 0) 



with 



dx 



(4) 



(5) 



I±ix)^TD. 



-1 ±V{x)/kBT 



,TV(y)/knT^y^ (6) 



Here xq is an arbitrary point and Dq— ksT/rj is the 
diffusion coefficient for Vq{x)— 0, i.e., for free diffusion. 
The particle current is given by [iJi 



(x) = {1 - e-^^/^--^) 



(7) 



The relationship between the directed and diffusive 
components in Brownian motion (the level of coherence) 
can be characterized by the Peclet number Pe or by the 
factor of randomness Q 



Pe 



D 



(8) 



The greater the Peclet number, the greater the coherence 
of Brownian transport. 



We proceed form the exact analytical formulas for the 
effective diffusion coefficient (0J and current iQ of Brow- 
nian particles, and consider a piecewise linear periodic 
potential with the amplitude A, period L, and asymme- 
try parameter k {0<k<L; the potential is symmetric if 
k=L/2). If the tilting force exceeds the critical value 



F>Fr,= 



A 



L^k' 



(9) 



the potential is monotonically decreasing, i.e., the poten- 
tial does not have local minima. 

We can take with no loss of generality L = l and re- 
place the relevant quantities with the corresponding di- 
mensionless ones: T = kBTA~^ and F — F/Fc, hence 
Fc = l, D = Dt]A-\ Do = DovA-^ so that Do = f and 
(x) —rjA^^ (x). For brevity, in what follows we will omit 
the tilde signs above the symbols. 

Now the expressions for the diffusion coefficient, cur- 
rent, and Peclet factor have the following form: 



D = TYZ-^ . 



Pe = ^oZ^iTY)-' 



where 



(^0 = 1- exp 



F 



T{1 -k)J ' 



Z= ( H_a{x)dx+ [ H^b{x)dx, 

Jo Jk 

y= / H+a{x)[H^a{x)]^dx 
JO 

+ [ H+bix)[H^b{x)fdx. 

Jk 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



Equations H14|) and (|15l) contain the functions H±a and 
H±b where the subscripts a and b associate, correspond- 
ingly, with the limits of integration from to fc and from 
fc to 1. We have 
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x+l 



^TVa{y)/T 



X k 
1 fc+1 



1 

x+l 



(16) 



k+l 



In Eqs. p6|l the dimensionless potential reads as 

Va (x) = flg — ax , < x< k, 
vi,{x) ^ —bQ + bx , k< x< 1 , 

Va (x) — ao — ax , 1 < x< 1 + fc, 
Vi^{x) = —bo + bx , 1 + k< x< 2, 



(17) 



with 



1 - {l~F)k 
(1 - k)k 

b = , bo = 



ao = 1 , do 



1 + k 



k 



bo 



1 + k 



1 - fc ' " 1 - fc ' 1 - fc 

Performing integration in Eqs. pti|l . we obtain 



(18) 



H±a{x) 



H±b{x) 



9fa 



a 

X exp 
— ^ + 



i[t2x - fc(l =F 1)] 
2T 



X exp 



with the notations 



6[±2(x-l) + (l-fc)(l±l)] 

2r 



1 1 



9 = 



1 - F 



b ' 



(19) 



(20) 



^Pa = exp 



T 



(^6 = 1- exp 



1 - (1 -F)k 



(21) 



1 - k 



:Tg'^fo 



1 1 



gifoVa^b 



62 



ffc 1 — fc . 

+2.gVo -<y3a(l-¥'&) ; — fbi'^ + Va) 

\_a b 

+Tg^ ^^liPb{l - Pb) + ^^IM^ + Va) 



(23) 



where 



(pa = exp 



2(1 -F) 



T 



"^6 = 1- exp 



2[1-(1-F)fci \ 

r(i-fc) ; 



(24) 



By that we have derived the exact algebraic expressions 
for the current (i), the diffusion coefficient D, and the 
Peclet factor Pe. 



IV. ASYMPTOTIC LIMITS AND PARTICULAR 
CASES 

In this Section we will examine the asymptotic limits 
and essential particular cases on the basis of the analyt- 
ical formulas derived. 

(i) In the absence of tilt (i^ = 0), Eqs. lO and if^ 
reduce to 

^0=0, (^, = £1/^-1, pb^l-e-'^^, (25) 
and from Eqs. lfTT|l . (|^ . and if^ one obtains 

1 



D 



2T [cosh(l/T) - 1] 



(26) 



This expression is as a special case of the general for- 
mula of the diffusion coefficient for arbitrary unbiased 
Substituting the functions H±aAx) from Eqs. ^ into periodic potential 23]. It is to be noticed, that for F = 



(|14|1 and p5|) . we have after integration 



a b 



ifO + Tg (paifib : 



(22) 



the coefficient of diffusion becomes independent of the 
asymmetry parameter fc. 

(ii) In the high temperature limit, one can take into 
account only the first order terms in the expansions of 
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the exponents in Eqs. (|13|l and H21() . Then 



F l~F 



r(i - k) 



T 



1 - {1-F)k 

m-k) 



and 



H±a.,b = T'\ Z = T-\ Y = T- 



(27) 



(28) 



The current, diffusion coefficient, and Peclet factor now 
become 



(i) 



F 



1-fc' 



D^T, 



Pe 



T(l-fc) 



(29) 



(30) 



(31) 



(iii) Under the conditions F > 1 and F/T > 1, it is 
valid that 



^0 ~ -•-Pa 1 

and 

i^±a.b 



1 - k 



1 - k 



F 



F 



-b: 



Y- 



F 



1 - k 



1 - k 



F 



(32) 



(33) 



As a result, the expressions for £), (i), and Pe coincide 
with Eqs. (|29|I - H31|I . Thus, at high temperatures and at 
large values of tilting force, the transport properties of 
Brownian particles are the same. 

(iv) If 1 and (1 - F)/T>1, we have the following 
asymptotic limits: 



Then Eqs. CH), IIU, and lEl) yield 

r=iTgV(i-^)/^[(a-i-fe-i) 
+2(a-ie-^/(i-^)^ + 6-i)] 



•fa 



and 



D = 



27^g2g(l-F)/T ' 



(34) 
(35) 

(36) 
(37) 

(38) 



Pe = 



2ipo 



— 2tanh 



^0 



2T{1 - k) 



(39) 



(cf. also 0). If, additionally, the condition F/T(l-k'p>l 
is fulfilled, it is valid that e~^^"^''^~'^'«0 and ipo~l. Conse- 
quently, in the present case we have 2D= (x) and Pe—2. 
This indicates that an extremely exact stabilization of 
the level of coherence of Brownian transport occurs in 
this region of parameters. 

(v) At the critical tilt (F=l), it is valid that 



a I 

H±b{x) = — + — - 
a 1 



k /a[=F2a;-A:(l=Fl)l 
-^"P ^ 



cxp I - — 



y^Q[±2(a;-l) + (l-fc)(l±l)] 
2T 



(40) 



whereas in the low-temperature limit Eqs. (|10|l - p5|l and 
(PI yield 



2T 



{l-kf 



D 



2T 



i{i-ky 



Pe ■ 



(41) 



(42) 



(43) 



We observe that, for F— 1 the Peclet factor is constant 
and depends neither on the temperature nor on the asym- 
metry parameter. 



V. NON-MONOTONIC BEHAVIOR OF 
DIFFUSION AND COHERENCE 

The expression of the diffusion coefficient as a function 
of F given by Eqs. ((TTl), 1221), and reveals a quali- 
tatively similar behavior to that found in Refs. [12,13], 
exhibiting a resonant-like maximum if the temperature is 
sufficiently low. This effect is strongly influenced by the 
shape of the periodic potential, as illustrated in Fig. 1 
This plot shows that for positive bias (F > 0), the increase 
of the value of k favors the amplification of diffusion com- 
pared to free thermal diffusion. 

Next we consider the dependence of diffusion on tem- 
perature. Figure 2 displays the typical behavior of the 
diffusion coefficient as a function of temperature and 
asymmetry parameter for F= 0.95. The non-monotonic 
behavior occurs if the value of tilting force lies within a 
certain range Fi<F <F2 where Fi^2 depend on k. One 
can see that the maximum of D(T) becomes rapidly nar- 
rower and higher as k approaches unity. There exist the 
limiting values ks w 0.8285 and Fq « 1.1292: for k<kE 
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FIG. 1: The plot of the difTusion coefhcient vs the tilting force 
and the potential asymmetry parameter at fixed T — 0.1. 



5z 




FIG. 2: The plot of the diffusion coefficient vs the tempera- 
ture and the potential asymmetry parameter at fixed i^ = 0.95. 



or F> Fq the diffusion coefficient is a monotonic function 
of temperature for all values of bias and asymmetry co- 
efficient. Therefore, tfie {k, F)-space is divided into two 
domains where the analytical properties of the diffusion 
coefficient as a function of temperature are qualitatively 
different. 

The curves of the Peclet factor versus temperature are 
depicted in Fig. 3. The function Pe{T) passes through 
a maximum (curves 1-5) for F<Fc, which is also present 
slightly above the critical tilt (curve 6). With a further 
increase of F, the maximum of Pe(T) disappears. As 
seen in Fig. 3, the optimal level of Brownian transport, 
determined by the maximal value of the Peclet number, 
is sensitive to the shape of periodic potential: at F<Fc 
the optimal level of transport rises with an increase in 
k. At the same time, the larger values of k make a min- 
imum of D(T), which follows a maximum of D{T) at a 
higher value of temperature, deeper (the effect can be 
anticipated in Fig. 2). In this sense the situation is anal- 
ogous to the results of Ref . j3| i where the enhancement 
of the coherence of Brownian motion in a certain region 
of temperature due to frictional inhomogeneity is associ- 
ated with the suppression of diffusion by the same factor. 

Now we discuss the amplification of diffusion by bias in 
terms of probability distribution. The stationary prob- 
ability density for the coordinate of a Brownian particle 




log(T) 



FIG. 3: The plot of the Peclet factor logj^Q(Pe) versus temper- 
ature log]^Q(r) for various values of the potential asymmetry 
coefficient at _F = 0.95 (curves 1-5) and F — 1.05 (curve 6). 
Curve 1: fc = 0.1, curve 2: fc = 0.5, curve 3: fc = 0.9, curve 4: 
fc = 0.95, curves 5, 6; fc = 0.99. 
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FIG. 4: The plot of the probability density P{x) at r = 0.01, 
k = 0.5 and F = 0.8. 



[l| normalized over one period can be written as 

^ J X 

Figures 4-6 represent (in terms of the dimensionless pa- 
rameters) the probability distributions characteristic to 
various diffusion levels depending on the tilting force. 

Figure 4 illustrates the situation, where particles are 
mainly localized around the minima of the potential and 
transport is strongly suppressed. Diffusion is essentially 
weaker in comparison with free diffusion: D/T r^lO~^. 
The probability distribution shown in Fig. 5 corresponds 
to the case where the diffusion is approximately maxi- 
mal {D/T— 3.5) for the chosen values of temperature 
and asymmetry coefficient. In this case the regions with 
a large probability are separated by the domains where 
the probability is much smaller, however, large enough 
to allow the entrance of a sufficient number of particles 
into these domains. As a result, a channel of hopping-like 



0.5- 



0.2 0.4 0.6 0.8 



1.2 



0.5 



1 1.5 

X 



2.5 3 



FIG. 5: The plot of the probabihty density P(x) at r = 0.01, 
k = 0.5 and i^=l.l. 



FIG. 7: The comparison of the dependencies of the Peclet 
factor and diffusion coefficient on the tilting force for various 
temperatures at fixed A; = 0.5. Sohd lines: 5xPe vs F, dashed 
lines: \o^iq[D{F) / D{Q)\ vs F. Curves 1: T^O.Ol, curves 2: 
r = 0.03, curves 3: T = 0.09. 



1 
P 

0.5- 




0.5 1 1.5 2 2.5 3 

X 

FIG. 6: The plot of the probability density P{x) at r = 0.01, 
fc = 0.5 and F = 3. 



transport is formed, leading to the enhancement of diffu- 
sion with respect to free diffusion. The further increase 
of the tilting force F makes the probability distribution 
still more homogeneous, as seen in Fig. 6, and the diffu- 
sion approaches the free diffusion limit {D/T~1.3 for the 
values of the parameters used in Fig. 6). 

Consequently, the amplified diffusion in the tilted peri- 
odic potential is characterized by the specific inhomoge- 
neous probability distribution with spatially alternating 
domains of high and low probability. The occurrence of 
a maximum in the dependence D(T) can be understood 
in a similar way. 

VI. CORRELATION BETWEEN THE 
ENHANCEMENT OF DIFFUSION AND 
STABILIZATION OF COHERENCE 

With regard to the simultaneous enhancement of diffu- 
sion and current, caused by the force F, with respect to 
an untilted system, the relation between D and (x) is of 
interest. One can expect that such a relationship reflects 
some intrinsic features of the mutual influence between 



FIG. 8: The comparison of the dependencies of the Peclet 
factor and diffusion coefficient on the tilting force for fc = 0.1 
(solid lines) and k = 0.9 (dashed lines) at fixed T = 0.01. 
Curves 1: logio[D(F)/D(0)] vs F, curves 2: 5xPe vs F. 



diffusion and current driven by the tilt merely at lower 
temperatures, when the initial suppression of both com- 
ponents of Brownian transport by periodic potential is 
stronger. 

The comparative plot of D and Pe versus F is pre- 
sented in Figs. 7 and 8. One can see that the function 
Pe{F) has a point of inflection which turns into a wide 
plateau at lower temperatures. For values of F, from zero 
up to the end of the plateau, the behavior of Pe{F) is de- 
scribed with great accuracy by Eq. H39|l . Note especially 
that the domain where Pe = 2 coincides with the domain 
where the increase of diffusion coefficient as a function of 
F is most rapid and follows with a good approximation 
the law D ~ c°'i+°'2F ^ where c and ai_2 depend on T and 
k. The location of the end of this region at larger values 
of F is quite insensitive to the shape of the periodic po- 
tential, as seen in Fig. 8, and is located approximately at 
Fc. 
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VII. CONCLUSION 

In the present work we have addressed the problem of 
overdamped motion of Brownian particles on tilted piece- 
wise linear potentials in the presence of white thermal 
noise. We have derived the exact algebraic expressions 
for the diffusion coefficient, particle current, and factor of 
randomness, which are valid for an arbitrary simple saw- 
tooth potential, tilting force, and strength of the thermal 
noise. On the basis of these formulas we have discussed 
the characteristic limiting cases and asymptotic behav- 
ior. 

Acceleration of diffusion turns out to be very sensi- 
tive to the shape of the piecewise linear potential. It is 
shown that large values of the asymmetry parameter k 
favor the amplification of diffusion by means of biased 
potential and temperature in comparison with free ther- 
mal diffusion. This can be understood as a result of the 
formation of probability distribution with spatially al- 
ternating regions of specifically balanced high and low 
probability. The necessary and sufficient conditions for 
the non-monotonic behavior of the diffusion coefficient as 
a function of temperature are established. On the basis 
of these results, among the other applications, one can 
expect the significant increase of diffusive transport, for 
example, in the superionic conductors p3.[25l| in a strong 
external electric field. 

We have demonstrated that the Peclet factor as a func- 
tion of the tilting force has a plateau at low temperatures 



which terminates with a steep rise at the critical value 
of the tilting force. As temperature grows, the plateau 
gradually reduces until disappears and the Peclet factor 
becomes monotonically increasing. The domain, where 
the Peclet factor exhibits the plateau, coincides with the 
domain where the enhancement of diffusion coefficient 
is maximal. Consequently, in the region of parameters 
where substantial acceleration of diffusion occurs, the 
current and diffusion are exactly synchronized. Note that 
the stabilization of the coherence level of Brownian trans- 
port takes place at the value of the Peclet factor Pe — 2 
which is also characteristic for Poisson processes [ p , su ch 
as the Poisson enzymes in kinesin kinetics l27l E^l . 

Finally, it is remarkable that, as a result of the inter- 
play of periodic potential, bias and white noise, there ex- 
ists an exact correlation between the acceleration of dif- 
fusion, induced by the tilting force, and the stabilization 
of the coherence level of Brownian motion. It seems that 
this phenomenon is quite universal and manifests itself 
for an arbitrary periodic potential in situations where ini- 
tially strongly suppressed transport is enhanced by bias, 
which generates significant amplification of diffusion in 
comparison with free diffusion. 
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